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C$ . Abstract 

In this paper, following Nachman's idea |Nac88j and Haberman and Tataru's idea 
[HTj . we reconstruct C 1 conductivity 7 or Lipchitz conductivity 7 with small enough 
C^l 1 value of |VZoi?7| in a Lipschitz domain Q from the Dirichlet-to-Neumann map A 7 . In 

the appendix the authors and R. M. Brown recover the gradient of a C7 1 -conductivity 
1 at the boundary of a Lipschitz domain from the Dirichlet-to-Neumann map A 7 . 
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1 Introduction 

Let C K n , n > 2 be an open, bounded domain and let 7 be a strictly positive real valued 
function defined on Q which gives the conductivity at a given point. Given a voltage 
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potential / on the boundary, the equation for the potential in the interior, under the 
assumption of no sinks or sources of current in 0,, is 

(1.1) div(7Vn) = 0, in O, u\ an = f. 
The Dirichlet-to-Neumann map is defined in this CclSG clS follows: 

A 7 (/)=7§^n, 

where is the outward normal derivative at the boundary. For 7 G L°°(Q) and being 
the boundary of O Lipschitz, then A 7 is a well defined map from H^{dVt) to H~z(d$l). 

The Calderon problem concerns the inversion of the map 7 — > A~, i.e., whether A 7 
determines 7 uniquely and in that case how to reconstruct 7 from A 7 . 

For the boundary determination: if dtt is C°° and 7 G C°°(Q), Kohn and Vogelius 
|KV84| showed that A 7 determines for all k > 0; if <9$7 is Lipschitz and 7 G Lip(Q), 

Brown |Bro01| showed that 7 Ian can be recovered from the knowledge of A 7 ; if dQ is C 1 
and 7 G C 1 (fi), Nakamura and Tanuma [NTOlj showed that g^|sn can be recovered from 
A~. Now, being 7 G C 1 (fi) and for Q Lipschitz domain, the authors and Brown together 
show the recovering of the gradient of the conductivity at the boundary (see Appendix). 

For the higher dimensional problem (n > 3) , uniqueness was first proven for piecewise 
analytic conductivities by Kohn and Vogelius |KV85j . Later, Sylvester and Uhlmann 
showed that the uniqueness holds for smooth conductivities in their fundamental paper 
[SU87J which opened the door for studying the Calderon problem. Generalization to 
less regular conductivities had been obtained by several authors. Brown |Bro96j obtained 
uniqueness under the assumption of |+e derivatives. Paivarinta, Panchenko and Uhlmann 
[PPU03| proved uniqueness under the assumption of | bounded derivatives. Brown and 
Torres |BT03j obtained uniqueness under the assumption of | derivatives being in L p ,p > 
2n. Later, Uhlmann jl.'lilOi) proposed a conjecture whether uniqueness holds in dimension 
n > 3 for Lipschitz or less regular conductivities. 

Recently, Haberman and Tataru jHTj gave a partial answer to this conjecture. They 
used a totally new idea to construct CGO solutions in Bourgain's space and showed unique- 
ness for C 1 conductivity 7 or Lipchitz conductivity 7 with small enough value of |VZo<77|. 
The ideas and techniques coming from this work have been widely used as can be seen in 
the papers by Zhang [Zhal2| . Caro, Garcia and Reyes [CGR12j or Caro and Zhou [CZ12] , 

If uniqueness holds, one whether or not construct the conductivity in the domain 17 
by the boundary measurements. Nachman |Nac88j and Novikov |Nov88] independently 
solved this problem for C 2 conductivity. 

We will briefly describe Nachman's idea as follows. For 7 G C 2 (Cl), if u is a solution 
to (11. ip . Sylvester and Uhlmann reduced v = -y^u to satisfy 

(1.2) -Av + qv = in ft, 

where q = ^l/ 2 , an d the corresponding Dirichlet-to-Neumann map A q to (II. 2j) is deter- 
mined by A~. For • p™ = 0, p^f 1 G C n and \p < if > \ — > 00, as n — > 00, Sylvester and 
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Uhlmann constructed a family of CGO solutions = e x ' p i (1 + <&^(x)) to 
(1.3) -Av + qv = inR n , 

with the remainder term decaying to zero in some sense as \pi\ — > oo, where 

q = ^1/2 in ^ and q = outside f2. 

For the appropriate chosen ' P2 = 0> /°2 ^ ^ n an d i°i + P2 = Green's 
formula gives us 

(n) 

(1.4) 



(x)e ia>fc (l + $£ n) (a;))ds = / (Age 3 "^" 5 - ge ^ 2 )v in) dS, 

Jan ou 



where 



(1-5) A ? =( 7-^-1/2 + i 7 - 1 ?J 



2 ' <9z/. 

Letting n — )• 00 in (ll.4p we conclude from the decay property of 

(n) 

/" (n) 5e X ' P 2 , \ 

(1.6) q(k) = lim / (A q e x ^ - — )v {n) dS, 

n ->°° Jan (Jv 

so the problem is then to recover the boundary values v^\qq = /(«). 

Pi 

On the othe hand, v ^ is a solution to the exterior problem: 
f -A«W =0 in R n \H, 

^• 7 ) 1 „,(n)| _ f 9»W| _ A f 

[ « v J \an - J p (n), —fcr\an — A ?/ p (™)- 
If w" satisfies the radiation condition 

(1.8) lim / G i„)(x,y) \ * - g („(") - M )dS(y) = 0, 

using Green's formula in (|1.7|) we can reduce / („) to satisfy a boundary integral equation 

Pi 

(1.9) /.w = e*'"? - (S (»,A, - B (n) - J/)/ ( „, 

Pi Pi Pi Z Pi 



and (jl.9p is uniquely solvable by Fredholm alternative theory, where G („) , S („) and B („) 

Pi Pi Pi 

are defined in Section 3. 

Unfortunately, it is not easy to directly check that u (n) satisfy ([Tgj) . To move around 
this obstacle, Nachman |Nac88j dealt with this problem in an inverse direction: He first 
constructed CGO solutions to (II. 3D automatically satisfying the radial condition (|1.8p from 
boundary integral equation (jl.9p and then showed that these solutions are the same as 
the solutions obtained by Sylvester and Uhlmann [SU87] . 
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For 7 G C 1 (f2) or 7 G Lip(Q) with the small enough value of \Vlog^\ we will follow 
Nachman's idea to construct CGO solutions to the equation 

(1.10) div( 7 Vu) = 0, in R n 

from the boundary integral equation (jl.9p . In view of the less smooth regular 7, we need 
to do some changes in the above steps. First we reduce conductivity 7 to the case 7 = 1 
near the boundary of -Br(O) D O and show that the new Dirichlet-to-Neumann map A 7 
corresponding to 7 in Br(0) is determined by A 7 . Next we construct CGO solutions to 
(jl.lOP from the boundary integral equation (|1.9p on the boundary dBn(0) and show that 
these solutions are the same as the solutions obtained by Haberman and Tataru |HT| . 
We state the theorem as follows. 

Theorem 1.1. The conductivity 7 can be recovered in the domain O by the knowledge of 
A 7 under one of the following assumptions: 

(a) Let Q C M n ,ra > 3 be a bounded domain with Lipschitz boundary and let j(x) G C 1 (0) 
be a real valued function with 7(0;) > Co > 0. 

(b) Let O C K n , n > 3 be a bounded domain with Lipschitz boundary and let "f(x) G Lip(£l) 
be a real valued function with 7(3;) > Co > and such that there exists a constant e n ^ 
satisfying \Vlog~f{x)\ < e n< n- 

Our paper is organized as follows. In Section 2 we reduce the conductivity 7 to be 
7 = 1 near the boundary. In Section 3 we construct CGO solutions from the boundary 
integral equation. In Section 4 we state the reconstruction of the conductivity 7 from 
A~. In the appendix the authors and R. M. Brown give the proof of the recovering of the 
gradient of a ^-conductivity at the boundary of a Lipschitz domain. 

Acknowledgements. The authors want to express their gratitude to Professor Mikko 
Salo and to Pedro Caro for many useful suggestions and deep comments. The authors 
would like to sincerely thank Professor R. M. Brown for pointing out the result appearing 
in the appendix. The first author is supported by the projects ERC-2010 Advanced Grant, 
267700 InvProb and Academy of Finland (Decision number 250215, the Centre of Excel- 
lence in Inverse Problems) and also belongs to the project MTM2007-62186 Ministerio de 
Ciencia y Tecnologia de Espaha. The second author is supported by Academy of Finland 
(Project number 2100001796, the Finnish Centre of Excellence in Inverse Problems). 

2 Reduction to The Case 7 = 1 near The Boundary 

For a bounded domain Q with Lipschitz boundary and 7 G Lip(fi), if A 7 is known, we 
can recover the values of 7 at the boundary dQ (see |Bro01| ). From this knowledge we 
can extend 7 to be a Lipchitz function in M n with 7 = 1 outside a large ball Br and 
l{x) > Co > 0. Furthermore, If \Vlog r y(x)\ < e n ^ we can keep this property holding in 
W 1 . For a bounded Lipschitz domain Q and 7 G C 1 (Q), if A 7 is known, we can recover the 
values of 7 and the gradient of 7 at the boundary d£l (see the appendix). Next, Whitney's 
extension allows us from the information of d a 'y(x) for all x G d£l and all \a\ < 1 to extend 
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7 to be C 1 in M. n with 7 = 1 outside a large ball Br and j(x) > Co > 0(The readers are 
referred to see Chapter VI §4.7 of |Ste70| ). In the above two cases obviously 7 is known 
in R n \n. 

Now for fixed R > R$ such that f2 C .Br(O), we define the new Dirichlet-to-Neumann 
map as follows: 



A 7 : f\dB R {o) — > inr\dB R (p), 



where uj is a solution to 



I div( 7 Vn)=0, in5 fl (0), 
1 u \dB R (0) = /■ 

In the following, we will show A 7 is determined by A 7 and the value of 7 in Br(0)\Q. 
Here we need to mention that Nachman already used this idea in |Nac96| and obtained 
an exact formula in the plane with the the conductivity 7 G W 2,p ,p > 1 (see Proposition 
6.1 of [Nac96| ). Following Nachman's idea, we will generalize this formula in a bounded 
Lipchitz domain Q C M. n ,n > 2, for a Lipchitz conductivity 7. 

In fact, if we denote by Oi, respectively the domains 0, and Br(0), and by <9f2i, 
respectively 8^2, the boundary of £1 and the boundary of Br(0), the Dirichlet-to-Neumann 
map corresponding to 7 in the domain f^X^i can be viewed as 2 x 2 matrix of operators 
A ij ,H 1 {dn i ) H5(dnj),i,j = 1,2., defined as follows. Given ft G H l {d£lj) for j = 1,2., 
considering the Dirichlet problem 

f div^Vw) = 0, in fi 2 \fii, 

\ MdQi = fl, w\dn 2 = /2, 

we have the following Dirichlet-to Neumann mapriT (dClj) — > L 2 (dflj) as follows: 



(2.1) 



7|^lan 2 rVA 21 A 22 A/: 



'2 

where J^la^ denotes the nontangential limit of from outside of f2i and J^lar^ 
denotes the nontangential limit of ^ from inside of Since 7 is known in f2 2 \Oi, 
A 1 - 7 , i,j = 1,2., are determined. Now for any / G H l (d£l2) we have a solution u G 
Fi(ft 2 )nC s (^2), where < s < 2, satisfying 

{div(7Vn) = 0, in fi 2 , 

Let g = u\qq 1 G // 1 (9f]i). We have that (|2.ip implies 

(2.2) A 7 = 7^|an 2 =A 21 < 7 + A 22 /. 
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If g can be recovered from A 1 - 7 , A 7 and /, then A 7 is known. On the other hand, in view 
of u G H^(Q 2 )r\C s {^2), we can deduce from (pTTj) 

A 11 p + A 12 / = -7|^l^ 1 = -A 75 . 

Hence, it follows that 

(-A 7 -A n ) 5 = A 12 /. 

If we can show that (— A 7 — A 11 ) is an invertible operator: ^(dfli) — > L 2 (dQi), then we 
have 

(2.3) ff = (-A 7 -A 11 )- 1 A 12 /. 



Finally (J22D and ([ZSD imply 

(2.4) A 7 = A 21 (-A 7 - A 11 )~ 1 A 12 / + A 22 /. 

Now we state the generalized result of Nachman |Nac96j as follows. 

Theorem 2.1. Let f^i d 0,2 be a bounded Lipschitz domain in R n , n > 2 and let ■j(x) G 
Lip(Q, 2 ) with j(x) > C > 0. Then (-A 7 - A 11 ) is an invertible operator: ^(dQi) — > 
L 2 (dfti) and A 7 = A 21 (-A 7 - A 11 )"^ 12 + A 22 . 

Before proving Theorem 12.11 we first introduce some known results. Let G(x, y) be 
the Green function such that, for iG^, 



div(7VG(x,y)) = 5 X , in fi 2 , 
G(x,y)| 9 n 2 = 0. 



(2.5) 

Now we define the single layer potential and double layer potential as 
(2.6) Sf(x)= [ G(x,y)f(y)dS(y), x G fi 2 , 



(2-7) D/(x) = / ^lllf(y)dS(y), x G n 2 . 

For the single potential Sf(x) and the double layer potential Df(x), we collect the fol- 
lowing results from Mitrea and Taylor's paper [MT99] (see Proposition 1.6, Theorem 3.7, 
Proposition 3.8 and Proposition 8.2 of |MT99| ). 

Proposition 2.2. If is a bounded Lipschitz domain in M. n ,n > 2 and f G L 2 (dQ,\), 
then the functions Sf{x) and Df{x) have the following properties: 

(a) The following estimates hold 

\\Df(x)\\ LHdai) <C\\f\\ L2(dUl) , 
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\\Sf(x)\\ HHdUl) <C\\f\\ L , idUl) . 

(b) div( 7 vs/(x)) = 0,16 n 2 \dn 1 . 

(c) Sf{x) G H x (Sl 2 - Oi) and G iT^Oi). The boundary value Sf+(x){Sf-{x)) 
of Sf(x) from outside (respectively inside) Qi are identical as elements of ^(dfli) and 
agree with Sf(x)\ dni . 

(d) The (nontangential) limits dS JJ^ i ^jf^ ) as the boundary is approached from out- 
side (respectively inside) Q,± are given by the formula 

dSfix) 1 1 

t; — = 7 --—fix) + B* fix), almost everywhere x G dU\, 

where 

B*f(x)=p.v. [ ^Af{y) dS {y). 
J dni du[x) 



In particular 



dSf(x) _ dSf(x) = 1 
dv^ dv + j(x) 



(e) The nontangential limits D f + (x)(D f_(x)) of Df(x) as we approach the boundary 
from outside (respectively inside) S7i exist and satisfy 

1 1 

Df+ = _ —fix) + Bfix), almost everywhere x G dQ\, 

27(a) 



where 

dG(x,y) 
du(y) 



Bf(x) =p.v. I 



f{y)dS{y). 



Remark 2.3. Our case is a little bit different from the case of the paper |MT99j . In |MT99j 
they dealt with the operator L = - A + V{x)(V{x) G L°°(M),V > 0,not identical zero.) 
in a compact, connected smooth M with C 1 metric tensor. But if we just want to get 
the above Proposition 12.21 we can reproduce the steps of the proof of Theorem 3.7 and 
Proposition 3.8 to show that the above Proposition 12.21 still holds for the operator L = 
div^Vu) in a bounded domain in with Lipchitz conductivity 7. 

We know that the Dirichlet-to Neumann map A 7 is a bounded operator:./? 2 (dQ\) — > 
H~z (dQ,\). In fact, when fix is a bounded Lipschitz domain and 7(2;) G Lip(Qi), the 
Dirichlet-to Neumann map A 7 is a bounded operator :H 1 (d^li) — > L 2 (<9r2i)(the same 
argument for A 7 ,A iJ ,z,j = 1,2.). This result follows from the following Proposition 12.41 
and Proposition 12.41 can be deduced from Proposition 7.4 and Proposition 8.2 of the paper 
|MT99| . 
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Proposition 2.4. Suppose ft is a bounded domain in M. n ,n > 2 with Lipchitz boundary 
and 7(x) G LipiQ) with 7(3;) > Co > 0. Ze£ u(x) G H 1 ^) be a solution to 

{div^Vu) = 0, in f2, 
u|an=5G^ 1 (^)- 

IIt^IIl 2 ^) < L\\g\\m(da) 

and 

\\(y u Y\\L2(dn) < C\\g\\ H i(Q Q ), 
where (Vn)* is the nontangential maximal fuction o/Vu. 

With the above preliminary works in hand we are now in the position to show an 
identity about A 7 + An. 

Lemma 2.5. For any g G H 1 (dQ,i), we have the identity 

5(A 7 + An)s'|aQ 1 = glddi.- 
Proof. For any g G ^(dQi) we can find a solution w G H 1 (Q,2\tli) satisfying 

J div(7Vw) = 0, in Q 2 \Qu 
\ Hani = 9, w\dn 2 = 0. 

Since 7(2;) G Lip(fl 2 ) the local regularity of the uniformly elliptic equation implies 
that w G H? oc {n 2 \di) and Proposition ELl implies that (Vw)* G L 2 (d(Q 2 \&i))- Therefore 
we can still use Green's formula for x G f^fii, 



w(x) = I di\{^Vw)G{x,y)dy — j w div (jVG(x,y))dy 

(2.8) =/ G (x,y)^dS(y) - [ ?§^ lw dS(y) 

Jan 2 dv Jdn 2 ov{y) 

+ [ G(x,y)A 11 gdS(y)+ [ ^^l 7wd S(y). 
Jaui Jani ou{y) 



Since G(x,y) =0,1/ £ dQ 2 and w\qq 2 = 0, fj2.8[) implies 

(2.9) w{x)=f G(x,y)A u gdS(y)+ [ 

Janj Jan! ov{y) 

Let x —7- dVL\ in ()2.9[) and Proposition 12.21 implies 



^wdS(y). 



(2.10) g(x) = [ G(x,y)A n gdS(y) + \g(x) + [ d( *^\ gdS(y),a.e.,x G d^. 
Jan! 2 J 9Ql du{y) 
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On the other hand we can find a solution u to 



[ div^Vu) = 0, in fix, 
\ = g. 

Using Green's formula again we know for x G fix 



u(x) = / div('jVu)G(x, y)dy — I ud\v{^VG{x,y))dy 

(2 11) 



dG(x,y ) 

an-L K ~ ,a '"^"" Kaj J dQi dv{y) 



G(x,y)A 1 gdS(y) - I a ..^\ ludS(y). 



Let x — > d£li in (12.111) and Proposition 12.21 gives 

(2.12) g(x) = ! G(x, y)A l9 dS{y) + \g{x) - [ >ygdS(y), a.e.,x G dn v 

Finally (pTL0l) + (pTL2j) imply 

g(x)= [ G(x,y)(A n +A 7 )gdS(y), 
JdQi 

which means 

5(A 7 + An)g|eni = gldfii- 

At last we need to point out that the above proof is just a formal deduction due to the 
singularity of div(7VG(x, y)). To overcome this obstacle we can deal with the integrations 
in a domain in which we remove a small ball B e (x) centered at the singular point x. Finally 
letting e — > we get the above identity. □ 

Now we will show that A 7 + An is an invertible operator: ii^dfii) — > L 2 (<9fii). From 
this Theorem 12.11 follows. 

Lemma 2.6. A 7 + An is an invertible operator: H 1 (dQi) — > L 2 (dQ.\) and the invertible 
operator is Sflg^ for any f G L 2 (<9fii). 

Proof. We first show that 5 , |qq 1 is an isomorphism from L 2 ((9fix) to i? 1 (9f2i). In fact 
Lemma 12.51 implies is surjective. To show that 5'|an 1 is injective we assume that 

f mi G(x,y)f{y)dS{y) = on dfii for / G L 2 (dfii). By Proposition EjZ] we know Sf(x) 
is a solution to div(7Vu) = in fix and Sf_(x) = on the boundary d£l±. Hence the 
maximum principle of the uniformly elliptic equation implies Sf(x) = in Oi. On the 
other hand Sf(x) is also a solution to div^Vu) = in fi2\fii. In view of the choice 
G(x,y) = 0,x G (9^2, we know Sf(x)\gn 2 = and Proposition 12.21 implies Sf + (x) = 
on the boundary dili, so the maximum principle of the uniformly elliptic equation also 
implies Sf(x) = in f^V^i- 
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From the analysis above we know that Sf(x) = in By Proposition 12.21 we have 

j_ as f{x) dsf(x) 
1 {x) I{x) dv_ du + 

which implies / = on the boundary dVL\ . Thus S\gsi t is bijective and Proposition 
12.21 implies S\qq, 1 is a bounded operator:L 2 ((9^i) — > Finally Banach inverse 

mapping theorem gives us is an isomorphism which L 2 (dQ,\) to H l (d£l\). 

On the other hand, from Lemma I2.5f s identity S\gn 1 (A^ / + An) = /, bijection of 
S'lgQj implies A 7 + An is also bijective and from Proposition 12.41 we have that A 7 + An 
is a bounded operator: ^(dtli) — > L 2 ((90i). Banach inverse mapping theorem implies 
(A 7 + An) -1 = S\qqi- □ 



3 Construction of CGO Solutions from the Boundary Inte- 
gral Equation 

In this section we will follow Nachman's idea to construct Complex Geometrical Optic 
solutions to the equation div^Vu) = from the boundary integral equation. We begin 
with some basic properties of Green's function. For p G C n with p-p = 0, (—^ 2 + 2ip-^)~~ 1 
is easily checked to be locally integrable as a function of £ in W 1 , n > 3. Its inverse Fourier 
transform formally written as 

g p (x)=g(x,p) = —j _ e + 2ip .^ 

is a tempered distribution satisfying 

(A x . + 2p ■ V x )g(x, p) = 5 (x). 

The distribution G p {x) = — e x ' p g(x,p) is then a fundamental solution for the Laplace 
operator in M n : 

-AG p (x) = 5 (x). 

It follows that G p {x) differs from the Green's function of classical potential theory by a 
global harmonic function in R n : 

Gp{x) = G (x) + H p (x), AH p = 0inR n , 

where 

G (x) = g(x, 0) = \x\ 2 - n (u; n = L^L). 

(n - 2)uj n r (n/2) 

Note that G p (x) is a smooth function for x away from the origin and has the same singu- 
larity near x = as that of Gq(x). 

Using the family G p of Green's function defined above we now consider analogues of 
the classical single and double potentials. Let 

S p f(x)= [ Gp(x,y)f(y)dS(y), 
Jan 
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where G p (x,y) =: G p (x — y) and to begin with, / is continuous on d£l and x G M n \dQ. 
Define the boundary, or trace , double layer potential by 

B p f(x) = p.v. [ 9G /*\ y) f(y)dS(y)M x G Ml. 
Jan ov{y) 

Now we collect some properties about single and double layer potentials from Nachman's 
paper |Nac88j as follows (see Lemma 2.4 and Lemma 2.5 of [Nac88j ) . 

Proposition 3.1. Let Q, be a bounded domain in ]R n ,n > 3 with C 1 ' 1 boundary and 
suppose f G Hz(dQ). Then the function u = S p f have the following properties: 

(a) Au = in R n \dn. 

(b) u is in H 2 (Q) and in H 2 (Q e ) for 6 > 9q, where 8$ is a large enough number so 
that Q C {x : \x\ < #o}> and for any 9 > 6q, £l e = {x : x ^ 0, |x| < 9}. 

(c) The boundary values u + {u-) of u from outside (respectively inside) Q, are identical 

3 

as elements of H^{dVt) and agree with the trace of single layer potential S p f. 

(d) u satisfies the following decay properties: 

(3.1) Um / G p (x,y) 7 ^-u(y) dG /*\ y) dS(y)=0,a.e.,x G R". 

o^°oj\y\ = o P dv(y) dv{y) 

3 

(e) Let g G Hz(d£l) and v = D p g defined in M n \<9f2. Then the nontangential limits 
v+ (v- ) of v as we approach the boundary from outside ( respectively inside ) Vt exist and 
satisfy 

1 

2' 



v+ = _-g(x) + B p f{x) a.e.,x G d£l. 



With the above preliminary works in hand we first deduce a boundary integral equa- 
tion for the CGO solutions. 

From the analysis we performed in Section 2 we can assume Q = B R (0) and j(x) = 1 
near the boundary dB R (0) and for 7 G Lip(B R (0)) (or 7 G C 1 (B R (0))) define the Dirichlet- 
to-Neumann map as follows 

A 7 : / G H*(dB R (0)) ->• —L G H2(dB R (0)), 

ov 

where ut is the solution to 

f div( 7 Vn) = 0, mB R (0), 
[ u \dB R (0) = /■ 

Clearly A 7 is a bounded operator: H2(dB R (0)) -»• H*(dB R (0)). 
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We are now ready to establish a one to one correspondence between the solutions of 
the boundary integral equation and those of the following exterior problem: 



(i) Aip = in R n \B R (0), 

(3.2) 



(it) ip G H 2 (B R ,\B R (0)), iorR'>R, 



(in) ip(x,p) — e x ' p satisfies (pTT 

M U = A 7^ on dB R (0). 

We state Lemma 2.7 of Nachman's paper |Nac88j as follows. 



Proposition 3.2. (a) Suppose ip solves f (|3.2p (i) — (iv)). Then its trace on the boundary 

fp = ^+= *l>\dB R (0)) solves 

(3.3) f p = e x " - (5 P A 7 -B p - h)f p . 

3 

(b) Conversely, suppose f p G Hz(dB R (0))) solves (|3.3p . T/ien the function ip(x, p) defined 
for x in R n \B R (0) by 

(3.4) iP(x,p) = e*-r-(S p A 1 -D p )f p 

solves the exterior problem <\3.2\) (i), (ii), (Hi), (iv). Furthermore ip\dB R (o)) = fp- 

In the following we will use the boundary integral equation (|3,3p to construct CGO 
solutions to div(7Vn) = in M. n . 

Lemma 3.3. Suppose f p G H^(dB R (0))) solves (J33J). T/ien 
faj There exists a unique solution u G -f^ 2 oc (]R n ) to 

(3.5) div( 7 Vu) = 

in M n suc/i f/iai u = ip(x,p) in R n \B R (0). 

(b) Let v = j2u and v G (W 1 ) is a weak solution to the following Schrddinger 
equation 

-Av + qv = 0,in R n , 

A 7 V2 

where q = — jy^— . Furthermore the following identity holds 

(3.6) v(x) = e x P + [ V-y^ 2 -V{-^G p {x,y)v)dy, a.e.,x£R n . 

(3.7) $(x,p)=l g p (x,y)q(y)dy + I g p (x,y)q(y)$(y, p)dy 



(c) Let v(x) = e x 'P(l + $(x,p)). Then $(x,p) G H} (W 1 ) satisfies 



as a tempered distribution. 
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Proof. For f p G H2(dBn(0))) we can find a solution w G H 2 (Br(0)) such that 

J divfrVto) = in B fl (0), 

l w \dB R (0) = fp- 

Proposition EE2 implies w\ 9Br ( 0) = ip(x, p)\dB R (o) = fp an d ^7 = g^f* ■ Hence defining 

j w in Bfl(0), 

[ iP{x,p) inM n \B fl (0), 

we know it G (R n ) is a solution to 

div(7Vu) = in R n . 

The uniqueness follows from the maximum principle of the uniformly elliptic equation. 
Then (a) follows. 

Letting v = 72 u in view of 7 = 1 in M n \.BR (0)(.Ro < R), we can deduce from f)3.5j) 
that v G Hf oc (M. n ) f| H? oc (M. n \B Ro (0)) is a weak solution to 

(3.8) - Av + qv = 0,in M n , 

where q = — 77^-- 

7-V^ 

For any R > R and 5 > small enough, we can consider a Lipschitz function hs(t) 
given by 

1 t€ [0, i?']\[-R — + 5], 



^ i G [R — 5, R], 



and then define the function H$(x) = h$(\x\) G Lip(B R i (0)) , x G B R /(0). Clearly, we have 

(i) < 1 and fl<y(x) | SSji (o) = 0, 

(it) S pt|Vff 5 (^)| CT j= : 5 i?+5 (0)\ J B R _ (5 (0). 

Given g > we define the following regularized version of G p (x,y) by 

(3.9) G e p (x, y) = L— (\x - y\ 2 + e 2 ) 2 ^ + # p (x - y). 

p (n- 2)uj n 

Taking the test function H 5 (y)G 6 p (x,y) G H^(B R (0)) in ([33]), we have 

(3.10) / VvV(H s (y)G p (x,y))dy= [ V 7 1/2 • V(-^vH s (y)G £ p (x, y))dy. 
Jb r (o) h Jb r (o) 7 i/2 
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On the other hand since v € H 2 (B R , (0)\B R (0)) and 7 = 1 in B R > (0)\B R (0), Green's 
formula gives us 



(3.11) 



[ V^ 2 -V(—vH s (y)G p (x,y))dy = f -AvH s (y)G p (x, 

r 7 r b Tr ' 

/ VvV(H s (y)G p (x,y))dy- / ^-G p (x,y)dS(y), 
Jt r , JdB R ,(0)OU 



y)dy 



where T R i = B R , (0)\B R (0). 

Now from ([3TT0]) + (f37TT|) we have 



(3.12 



V ^ /2 ■V(-^ H vH s (y)G £ p (x,y))dy 
fl„/(o) T 1 



(3.13) 



/ • V(H s (y)G p (x, y))dy - [ ^G p (x, y)dS{y) =: I + II. 

Jb r ,(o) JdB Rl {0)OU 

Note that the left hand side of (|3.12p can be rewritten as follows 

W /2 ■y(^vG p (x,y))H 5 (y)dy 

r 1 

/ V7 1 / 2 • VH 5 (y)— E vG £ p (x,y)dy =: III + IV. 
Jb r ,(o) 7 1/2 

Letting 5 — > Lebesgue dominated convergence theorem implies 

(3.14) lim 111= [ V^ 2 -V(-^vG p (x,y))dy. 

On the other hand in view of spt\V Hs(x)\ C T$ =: B R+ s(0)\B R -$(0) and 7 = 1 in 
R n \B Ro (0), when 5 < R - R we have 

(3.15) IV = 0. 

Combining (|37L2j) . (|3TT3j) . (1311) and (ETTBl . we deduce that 
(3.16) 

lim / V 7 1/2 • V(-^vH s (y)G p (x, y))dy = f V 7 1/2 • V(-^vG p (x, y))dy. 
For I we have 

(3.17) / = / Vv VG £ p (x, y)H s (y)dy + f Vv VH 5 (y)G p (x, y)dy =: V + VI. 

Jb r/ (o) Jb rI (0) 

For V by letting 5 — > Lebesgue dominated convergence theorem implies 

(3.18) limV= [ Vv -VG £ Jx,y)dy. 

Jb rI (o) 
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In view of spt\VHg\ C T$ and v G Hf oc (JH n \BR Q ) integration by parts gives us 

VI = I VvVH s (y)G £ p (x,y)dy 
(3.19) =-/ AvH s (y)G £ p (x,y)dy- [ Vv ■ VG £ p (x,y)H s (y)dy 



IT S JT S 



JdB R+s (o) du JdB R _ s (o) du 

Since u G ^(R^S^) f) C^R^-BrJ and G^(x,j/) is smooth, letting 5 -> in (1X131) 
we have 

(3.20) lim VI = 0. 

<5->0 

Combining (I3TT2D . (IBTTD . (^T8l) and (^201) . we have 

(3.21) lim / Vv ■ V(H s (y)G £ p (x, y))dy = [ Vv ■ VG £ p (x, y)dy, 
and from (fBTHj) . ([3TTB]) and (|3T2Tj) we derive 



(3.22^ "^' (0) 



/ W /2 -V(-^vG £ p (x,y))dy 
^ Vu-VG£Oc,j/)dy- / 

(0) J dB (0) 



^S fl /(0) JdB„,(0)O l/ 

On the other hand, Green's formula gives us 

vAG £ p (x,y)dy 



B R ,(0) 



(3.23) 

From (13321) and IET251) . we deduce 
(3.24) 



/■ /■ dG e Jx,y) 



dv f dG £ Jx,y) 

^G £ p (x,y)dS(y)- / ~ 1 P T 4 
b rI (o) p JdB R ,(o)dv p JdB B ,(o) dv{y) 



vAG £ Jx,y) = G zUy) dS {y)- p\ \ ' vdS{y) 



+ I V^ 2 -V(^vG £ p (x,y))dy. 



Due to 7 = 1 in M. n \B R i(0) and v = 7 1//2 u, letting e — > in (|3.24p we have for almost 
every x G B R > (0) 

v ( x ) = f TrG p (x,y)dS(y) - [ 9Gp udS(y) 
(3.25) J8B Rl (0)dv J 9BrI{0) du(y) 

+ / V^ 2 -V(-!-vG p (x,y))dy, 

JR" 7 1 
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where we have used that —AG £ p (x,y) is an approximation of identity and L p estimate of 
the convolution type integration about the kernels of Gq(x — y) and VGq(x — y). 

In view of —A y G p (x,y) = 5 X and Ae x ' p = whenever p ■ p = 0, we deduce from 
Green's formula 

(3-26) e^=f ^l Gp (x,y)dS(y) - [ e xp dS{y). 
JdB R ,(o) ov JdB R ,(o) du{y) 

Since u = ip{x, p) in W 1 \Br{{)) and Proposition 13.21 implies that u — e x ' p satisfies (|3.ip . we 
can deduce from (13.26|) 

(3-27) lim ( [ ^G p (x, y)dS(y) - [ ^/*' y) udS(y)) = e xp . 

Then (13351) and (13371) imply 

(3.28) v(x) = e x - p + [ V 7 1/2 • V(-^G p (x, y)v)dy, a.e.,x G R n . 

Jr« 7 1 

Then (b) follows. 

Let 7 e = J e * 7, where J £ is the standard mollifier. Since 7 G LipiW 1 ) and 7 = 1 in 
M. n \Bn(0), by integration by parts we can deduce from f)3.28|) 

r a 1//2 

(3.29) v(x) = e x - p - lim / -^G p (x,y)vdy. 

Recalling v(x) = e xp (l + &(x,p)) we can deduce from (|3,29D 

f A. ^ 2 f A. ^ 2 

(3.30) $(z,p) = lim/ g p (x,y) — r^-dy + lim/ g p (x,y) — -^-<$>(y,p)dy, 

where we have used the relations G p (x,y) = G p (x — y), g p (x,y) = g p (x — y) and g p {x) = 
-e x -PG p {x). 

Hence, we conclude from (|3.30p that 

f A7 1 / 2 f A7V2 

(3.31) $(x,p) = g p (x,y) — - pr dy+ g P (x,y) — -j-<S>(y,p)dy 

Jw* 7 1/2 JR" 7 1/2 

as a tempered distribution. 
In fact, 

< 3 ' 32 > (1 «> = HlFT^n (7^) «>' 

and 

< 333 > (/„„ «•(*■ »)»<»)•<»■ (« - -|f I' + «/■•{ (l^ 5 *' «>• 

For 7 G Lip(R n ) and 7 = 1 in R n \B R (0) and $(x,p) G -ff^ oc (M n ) from the results of 
the paper of |HT| . we know (|3.31 j) . (|3.32p and H3.33P all make sense. Then (c) follows. □ 
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In Lemma llPl we proved the identity (|3,7p . We mentioned that J Rn g p (x, y)q(y)dy and 

Jwn 9p{ x -,y) ( l{y)^{y-: P)dy are tempered distributions. In fact in the paper |HT| by Haber- 

i 

mann and Tataru, they showed J Rn g p (x,y)q(y)dy G Xp 2 and J R „ g p (x,y)q(y) $(y,p)dy G 

i 

Xp . We are now in the position to introduce the definition of Bourgain's space X b p (see 
|Bou93] ) . Following the idea of Haberman and Tataru we define the space X h p by the norm 

where p p (£) = -|£| 2 + 2ip ■ £ is the symbol of A p := A + 2p ■ V. 
We can prove the following result. 

Theorem 3.4. Let 7(2;) G Lip(B R (0)) be a real valued function and assume that j(x) > 
Co > with j(x) = 1 near the boundary dB R (0) . Then there exists a constant e n)R such 
that if 7 satisfies either ||VZop7|| < e n ^ R or 7 G C 1 (5r(0)), then the following properties 
hold 

(a) S p A~f — Bp — is a compact operator on H2 (dB R (0)). 

f&j For any p ■ p = 0, p G C n , u>/ien |p| is /an/e enough, then there exists a unique 
f p G H^(dB R (0)) such that 

f p = e x ? - (5 P A 7 - B p - l -I)f p . 

(c) For any p- p = 0, p G C n , £/ien £/iere exists a sequence of p^ with \p( n '\ — > 00 such 
that = 7 - 1 /2 e ^-P ( " ) (i + $( Xj p(«))) are solutions to 

div(7Vu) = 0, in IP. 

Moreover 



|*(x,pW)|| 1 ->0, 

»(") 



.-i — > 0, 

X P (") 

as p*-" - ) — > 00, where t]b is a smooth function with compact support. 

Proof. For (a) let g G H^(dB R (0)) and consider it G H 2 (B R (0)) solution to 

f div( 7 Vn) =0, in 5 fl (0), 
i n l9B ii (0) = 5- 

Let x G and apply Green's formula to G £ p (x,y) in After taking the limit as e tends 
to zero we obtain 

/ G p (x,y)Audy + u(x) =[ G p (x,y)^dS(y) - [ G p (x,y) dGp ^ y) gdS(y) 

JB R (0) JdB R (0) OU JdB R (0) OV 

= {S p k 1 -D p )g{x). 
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Letting x —> dBfj(0) (nontangential from inside Br(0) ), we deduce from Proposition 13.11 

(3.34) (S P A 7 -B p - \l)g = T [ G p (x, y)Vlo gi • Vudy, 

2 Jb r (o) 

where T denotes the trace operator. 

Since 7 G Lip(B R (0)) from the estimates of the uniformly elliptic equation we know 
that the operator P 7 : g G H 2 (dBn(0)) — > u G H 2 (Br(0)) is bounded. Furthermore 
Rellich's compact embedding theorem implies that the operator S : it G H 2 (Bji(0)) — > 
V/o<77 • is compact and C alder on- Zygmund estimates imply that the operator G*p : 
/ G L 2 (Br(0) — >• G p / G H 2 (Br(0)) is bounded. On the other hand, the trace operator 
T : H 2 (B R (0)) -> H2(dB R (0)) is bounded. Summing up the above analysis we can deduce 
from (|3.34p that S p A 7 — -Bp — \l = TGpEP^ is a compact operator on H a ($£?#(())). Then 
(a) follows. 

To prove (6) by Fredholm alternative theorem we just need to show that the homoge- 
neous equation 

(3-35) /p = (_s p A 7 + 5p + i/)/p 

3 

only has the zero solution. For any g G H^{dBn{Q)) satisfying (]3.35p repeating the steps 
of the proof of Lemma [3.31 we can find a solution u(x, p) = i y~ 1 / 2 e x ' p &(x, p) G Hf oc (M. n ) to 

div^Vu) = 0, in R n , 

with u\ 9Br ( 0) = g and 

(3.36) $>(x,p)= g p (x,y)q(y)$(y,p)dy. 

Noting &(x, p) G Hj-^W 1 ) and q = A 1/ 2 with compact support, by dual method we can 
show 

(3-37) q{y)${y,p)exp. 

Then (1336]) and (|5^7|l imply $(x,p) G X} . 

Under our assumptions: 7 G Lip(B R (0)) and ||VZo<77|| < e n ^ or 7 G C x (Br(^S)\ when 
p is large enough Haberman and Tataru |HT| showed the operator : 

fix) G Xp > / g P {x,y)q{y)f(y)dy G X} 

is a contraction. Therefore it follows from (|3.36p that <J>(x, p) = 0. Hence u = and g = 0. 
Then (6) follows. 

In view of (6), for any p ■ p = 0,p G C", when \p\ is large enough, then there exists a 

3 

unique f p G H^{dB R {Q)) such that 

f p = e x ? - (5pA 7 -Bp- l -I)f p . 
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Now Lemma I3T31 implies that u(x) = 7 1 ' 2 e x ' p (l + <fr(x, p)) is a solution to 

div( 7 Vu) = 0, in R n , 

and 

(3-38) $(x,p)= g p (x,y)q(y)dy + g p (x,y)q(y)<S>(y, p)dy. 

Finally from (|3.38p under our assumptions about 7, Lemma 3.1 and Theorem 4.1 in |HT| 
imply that (c) follows. □ 

4 Reconstruction of the Conductivity 7 

Theorem 4.1. Under the assumptions of Theorem \3.4\ j can be recovered from the knowl- 
edge o/A 7 at the boundary 8Br(0). 

Before proving Theorem 14.11 we first recall an integral identity appearing in Lemma 
4.1 in |Knu06| . 

Proposition 4.2. Suppose 7 j S C 1 (fi),« = 1,2. andu\,U2 £ -ff 1 ^) satisfy V -(jiVui) = 
m O. Suppose further that u\ G satisfies V • ( 7 iV"Ui) = wrai/i ui = «2 on OO. 

TTien 

/• 1 1 1 1 /• 

/ (7i 2v 7| ~ 7| v 7i ) • V(uiu 2 )da; = / "jid v {u\ - u 2 )u 1 ds, 
Jn Jon 

where the integral on the boundary is understood in the sense of the dual pairing between 

H^(dU) and H~^(dSl). 

Remark 4.3. Proposition 14.21 also holds for 7 , G Lip(Q). 
Given k £ M n , we set 

Pi = s„77i + 1 I - + r n ??2 

(n) . . fk 

p\ = -s n r)i +1 - - r n n 2 

where 771,772 £ S" 1 " 1 satisfy (fe, 771) = (k,rj2) = (771,7/2) = and ^j- + r\ = s^. The vectors 

are chosen so that pf ■ p™ = 0, i = 1,2. and p± + p^ = i/c. According the choices of 
pf, i = 1,2. of |HT| we can take p(n) = p™ with \p™\ — > 00 in the choices of (c) of Theorem 
1521 

Proof of Theorem 14.11 

Proof. Taking the conductivities 71 = 7 and 72 = 1 and the solutions u\ = j~ l / z e x ' p i (1+ 

$>(x, p±)) and = e x ' p i in Proposition 14.21 we have 
(4.1) 



(n) 

/ -V7 1 / 2 -V(7" 1 / 2 e-^ ) e -^ ) (l+d.(x,p^))dx = / {A^-^-^^^dS. 

Jb r (0) JdB R (Q) OU 
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We have that u^\^ n ^ R ^ is a solution of the exterior problem (|3.2p such that 

In) i 
u\ IdBvto) — / («)• Since e x ' p i as R co), S ( n ), B („) are known and A 7 and 7 at dBn(0) 
rtv ' Pi rtv Pi ' Pi 

can be determined, Proposition 13.21 and (6) of Theorem 13.41 imply that the right hand side 

of (|4.ip is known. By the choices of + p^ = ifc we deduce from (|4.ip 

(4.2) 

(n) 

/ -V7 1 / 2 .V( 7 - 1 / 2 e"- fc (l + $(x,p^ ) ))dx= / (A 7 e^ n) - 7 a(e " P2 - )/ ( n )dS. 
Jb r (o) JdB R (0) OV Pi 

Let 7 e be as in Lemma 13.31 Now integration by parts in the left hand of (14. 2p gives us 

(4.3) lira / ^|V fe (l + Sfop^Jcfc = / (A 7 e^ n) ~ l ^T h f^dS. 
Clearly we have 

r a 7e 1/2 . h r A7V2 . , 

(4.4) lim / -^-e"* = / ~^^ k = q(k) 
^°Jb r (o) 7e 1/2 Jb r (o) 1 ' 

in the sense of the tempered distribution. 

Since || 7e — 7 ||h 1 (R") — > as e — > it holds 

(4.5) hm I ^le^^p^dx = (^jlj^tpB, <*>(*, pft 
£ ^°Jb r (o) 7e 1/2 T' 1 

1 . _i 

where ( ) means the dual pairing between X 2 (n) and X A and (ps € Cq°(W 1 ) with 

Pi Pi 

<Pb = 1 in B R (0). 

By (c) of Theorem 13.41 in Section 3 we deduce from (|4.5p that 

(4.6) lim lim / -^W^^x, p^)dx = 0. 



o^Uoo £ ^°Jb r (0) 7e 1/2 



Finally from ()4.3|) . ()4.4p and (j4.6|) we have that g(fc) is known as a tempered distribu- 
tion. By inverting the Fourier transform the potential q{x) can be recovered in W 1 . From 

the definition of q(x) = A J/ 2 we know w = log^/j G Hq(Br(0)) is a weak solution to 



(4.7) 



Aw + \Vw\ 2 = q, mB R (0), 
w\dB R (o) = 0. 



The nonlinear Dirichlet problem (|4.7p has a unique solution by the maximum principle of 
uniform elliptic equations and since we have already recovered q(x) we may construct 7 in 
Br(0) by solving the equation (|4.7|) . In other words 7 can be recovered by the knowledge 
of A 7 at the boundary dBji(0). Then Theorem 14.11 follows. 

□ 
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5 Appendix. 



Recovering the gradient of a C 1 -conductivity at the boundary 
C R. M. BrowrQ Andoni Garcia Guo Zhang 
c Department of Mathematics, University of Kentucky 
Lexington, Kentucky 40506, USA 

The goal of this appendix is to give a method for recovering the gradient of a coefficient 
from the Dirichlet to Neumann map. The argument is an extension of a method from 
earlier work of Brown [BroOl] . 

Throughout this appendix, we let SI be a Lipschitz domain as defined in |Ver84j . for 
example, and we let 7 denote a function on Cl that is continuous and satisfies the ellipticity 
condition for some A > 0, 

(5.1) A<7<A" 1 . 

We will use H s to denote the standard scale of L 2 Sobolev spaces. Given / £ H l l 2 (d£l), 
we may solve the Dirichlet problem 

J div 7VU = in 

1 u = f, on dQ 

We define 'ydu/dv, the co-normal derivative of u, as an element of H^^^dQ) by 

du f 
(7-^-» = / jVu-Vcfidy. 
ov Jn 

Here, (•,•): H~ l l 2 {d9.) x H 1 / 2 ^) ^ C is the bilinear pairing of duality and we use (ft to 
denote both a function in H l l 2 (d£l) and an extension into ft which lies in H 1 ^). Because 
u is a solution, the right-hand side depends only on the boundary values and not on the 
particular extension chosen for (p. 

We give a reformulation of a result of the author Brown for recovering the conductivity 
at the boundary |Bro01j , 

Theorem 5.1. Let Q be a Lipschitz domain. For almost every x G d£l, we may find a 
sequence of functions fx which are Lipschitz on dQ,, supported in {y : \x — y\ < N^ 1 / 2 }, 
satisfy 

ll/iv||^OQ) < CW S " 1/2 , < 8 < 1, 

and so that if 7 is a continuous function satisfying \5. and un is the solution of the 
Dirichlet problem for div 7V with data /at, and ip is a continuous function in £l, we have 

t/j(x) = lim / i[)(y)\Vu N (y)\ 2 dy. 

1 E-mail address: c Russell. Brown@uky.edu 
This work was partially supported by a grant from the Simons Foundation (#195075 to Russell Brown). 
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As an immediate consequence of the above theorem, we obtain recovery at the bound- 
ary and a stability result. The stability result for the recovery of the boundary values of 
a continuous coefficient in smooth domains was proved by Sylvester and Uhlmann [SU88, 
Theorem 0.2]. 

Theorem 5.2. Let 7 be continuous function in Cl and let x G d£l and Jn be as in Theorem 
\5.1\ then we have 

j(x) = lira (A 1 f N ,f N ). 

N— >oo 

As a consequence, we obtain the stability result: If 7 and 7 are continuous on Cl and 
elliptic, then 

117 - 7lU°°(n) ^ C|I A 7 ~ A 7ll£(// 1 /2(an),H- 1 /2(an))- 
Proof. This follows immediately from Theorem 15.11 □ 

If 7 is continuously differentiable in the closure of O, we have additional regularity 
of the solution ujy. Since the boundary values /at lie in we have that the full 

gradient of un lies in L 2 (dQ) and we obtain estimates for the non-tangential maximal 
function of the gradient (Vun)* 

(5-2) ll(V«Ar)*|| L2(9 n) < C\\f N \\ HHm) < CN 1 ' 2 . 

This result holds for elliptic operators in Lipschitz domains with C 1 coefficients and follows 
from the work of Mitrea and Taylor on manifolds with C 1 -metrics [MT991 Section 7]. Will 
we use the estimate (|5.2p to justify the use of the divergence theorem for expressions 
involving the gradient of a solution on the boundary. In our next result, we let VfU = 
Vu — denote the tangential component of the gradient of u. The following theorem 
gives a recipe for recovering V7 from the gradient of /jv on the boundary, the boundary 
values of 7, and the Dirichlet to Neumann map acting on /jv- 

Theorem 5.3. Let Q be a Lipschitz domain and let x and fx be as in Theorem \5.1[ Let 

a be a constant vector. We have 

a-V 7 (x)= lim / (7|V t / 7 v| 2 + -|A 7 / 7V | 2 )Q-^-2K((A 7 / 7 v)a- Vu N )da 

Proof. The rather mysterious expression on the right is obtain by rewriting the Rellich 
identity |Rel40j using -A 7 /at for the normal derivative. The following identity holds 
because the integrands are equal as functions 

[ (l\V t f N \ 2 + i|A 7 /jv| 2 )a • v - 2K(A 7 / JV a ■ Vu N ) da 
Jan 7 

= / 7|Viijv| 2 o; • v — 2^(7— a ■ V«at) da. 
Jan ov 
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Recalling that un is a solution, an application of the divergence theorem gives 



f du f 

/ 7|Vii7v| 2 a • v — 2^(7— a ■ Vun) da = / a-V7|V 
Jan du J n 



un\ 2 dy. 



We use the non-tangential maximal function estimate (|5.2|) to justify the application of 
the divergence theorem. With this, the theorem follows from the properties of um given 
in Theorem I5JJ □ 
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